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Abstract. Attosecond streaking is one of the most fundamental processes in attosecond 
science allowing for a mapping of temporal (i.e. phase) information on the energy domain. 
We show that on the single-particle level attosecond streaking time shifts contain spectral 
phase information associated with the Eisenbud-Wigner-Smith (EWS) time delay, provided 
the influence of the streaking infrared field is properly accounted for. While the streaking phase 
shifts for short-ranged potentials agree with the associated EWS delays, Coulomb potentials 
require special care. We show that the interaction between the outgoing electron and the 
combined Coulomb and IR laser fields lead to a streaking phase shift that can be described 
classically. 



1. Introduction 

The emerging field of attoscience [1~7] enables the investigation of electron dynamics as well 
as timing information of photoionization processes. Attosecond streaking has developed into 
a powerful tool to achieve temporal- resolution on the sub- 100 attosecond time scale [8-11]. It 
is based on a pump-probe setting with an extreme ultraviolet (XUV) pulse of a few hundred 
attoseconds duration serving as pump and a phase-controlled few-cycle infrared (IR) pulse as 
probe. Temporal information about the photoionization process can thus be mapped onto 
the energy axis in analogy to conventional streaking. First proof-of-principle implementations 
enabled the direct measurement of the life time of the Kv(3d~^) hole by Auger decay [8], the first 
measurement of the time-dependent electric field of an IR light wave [9] ; as well as measurements 
of delayed photoemission from core levels relative to conduction band states of a tungsten surface 
[10] and from different shells of atomic neon [11]. In the latter case, the measured streaking 
time shift is one order of magnitude shorter than the XUV "pump" pulse and two orders shorter 
than the oscillation period T of the IR probe pulse. 

The challenge in interpreting the obtained time shifts lies in disentangling the modifications 
of the observed streaking time shifts by the IR probe field from the intrinsic time shifts one is 
interested in. These intrinsic time delays are connected to the spectral variation of the scattering 
phases by the Eisenbud-Wigner-Smith (EWS) time delay operator [12, 13]. In this contribution 
we will apply the EWS time delay formalism to photoionization (half-scattering) and show with 



the help of numerical simulations that the EWS time delay and thus scattering phases become 
accessible by attosecond streaking as long as the potential is short-ranged. In particular, we 
study attosecond streaking of the release time of electrons in atomic photoemission [11, 14] by 
solving the time-dependent Schrodinger equation for effective one-electron systems. While the 
EWS time shift only contains the scattering phase from the single-photon dipole transition, 
in streaking setups an additional (non-perturbative, in contrast to RABITT, cf. [15]) IR field 
is present. To disentangle effects that stem from distortions of the initial state by the IR 
field we employ a "restricted ionization model" [16, 17]. In addition, we perform classical- 
trajectory Monte Carlo simulations to account for trajectory effects on the time shift resulting 
from the interaction between the outgoing electron and the combined Coulomb and IR laser 
fields. When both distortion effects are accounted for, the EWS time shift agrees well with the 
associated streaking time shift for short-ranged potentials which we will show for the example 
of the centrifugal potential. Unless otherwise stated, atomic units are used. 

2. Computational Methods 

The Schrodinger equation for a one-electron system in the presence of an external IR and XUV 
field is given by 

= H\ij) , H = H, + H,^ + H^^y, (1) 

where Ha is the atomic Hamiltonian with an effective one-electron potential, H^w = ^■-^xuv(i) 
describes the interaction with the the electric field Fxuv of the attosecond pump pulse and 
Hi^ = r • Fiji{t) the interaction with the IR probe pulse. Typical energies and intensities of the 
electric fields are fkoxw ~ 100 eV and /xuv < lO^^W/cm^ for the attosecond pump pulse, and 
/iWiR K 1.5 eV and /ir ^ 10^° - lO^^ w/cm^ for the probe field. 

Our computational method for solving Eq. 1 is based on the well-established pseudo-spectral 
split-operator method as described in [18]. We also employ a constrained Schrodinger equation, 
in which the interaction of the bound state with the probe field is suppressed. This so-called 
"restricted ionization model" (RIM) [16, 17] is based on the decomposition of the wave function 
tp{r,t) into an initial state part ^^{r) and a continuum part ^{r,t) 

\^|;{t))=e-'''^'\^o)+e'^'\m), (2) 

where £'o^'o(^ = -f/a(^*I*o(r) E = Eq -t-Wxuv- Inserting this ansatz into the Schrodinger 
equation (1), leads to the following (exact) expression describing the time evolution of the 
continuum wave function $(r, i), 

i-|l>(t)) = (h^ + Fir, + ^xuv - e) |$(t)) + [e'^^^^^'H,^, + H+^^ + e^^-xuvt^-^^^ j^^)^ (3) 

where -ffxuv = e*^^uviff-^^_|_g-f^xuvtjif+^^_ Following Schafer et al. [16], we apply the following 
simplifications: We neglect the interaction of the XUV pulse with the continuum wave packet 
and drop the highly oscillatory terms e^^'^^^vt^-^^ e^^^^^^Hi^ acting on the initial state 
(slowly varying envelope approximation) . The latter simplification suppresses the interaction of 
the IR pulse with the initial state, thereby inhibiting ground-state polarization and excitation. 
Their influence on the streaking time shifts can thus be investigated by comparing RIM and full 
TDSE results. The equation governing the time evolution of the continuum wave function then 
has the structure of a Schrodinger equation with an additional source term and can hence be 
solved in the same way as Eq. 1. 
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Figure 1: (a) Streaking spectrograms for an 800 nm IR laser field with a duration of 3fs and an 
intensity of lO^^W/cm^. The graphs show the final momentum distribution along the laser 
polarization axis for a H(ls) and He'''(ls) initial state and XUV photon energies of 80 eV 
together with the vector potential Ai^(t). The solid white lines are the first moments of the 
electron spectra. The dashed white vertical lines indicate the shift of the central minimum 
relative to the vector potential (red solid line), (b) Temporal shifts ts extracted from quantum 
mechanical streaking simulations with the full TDSE (full circles) and RIM (open diamonds), 
classical (CTMC) streaking simulations (open squares), and for comparison, the EWS time shift 
^EWS = dip/dE applied to the Coulomb phase (H"*": red solid line, He^"*": blue dashed line). 



3. Results 

We will start our discussion on the accessibility of scattering phases by attosecond streaking 
with an investigation of streaking time shifts for atomic photoionization in hydrogenic systems. 
The attosecond streaking technique is based on the assumption that the asymptotic momentum 
of an electron ionized by an ultrashort (few hundred attoseconds) XUV pulse is shifted by the 
instantaneous value of the vector potential j4ij^(r) of the IR streaking field at the moment of 
ionization, 

Pf{T)=PQ-A,^{T). (4) 

Here, po = [2(wxuv — ^i)]^ is the unperturbed asymptotic photoelectron momentum, where 
Wxuv is the XUV photon energy and £i the initial binding energy (all in a.u.). For temporally 
well-defined IR electric streaking fields Fi^{t) and vector potentials j4iR(r) = — FiYi{t)dt 
the momentum shift (energy shift) of the photoelectron allows thus, in principle, a mapping of 
emission time r onto energy ('streaking'). 

Fig. la shows photoionization spectrograms for ionization from the hydrogen and He^ ground 
state by an 80 eV XUV pulse with a duration of 200 as (FWHM of the Gaussian intensity 
envelope), streaked by an 800 nm IR laser field with a duration of 3fs (FWHM of the cos^ 
envelope) and an intensity of lO^^W/cm^. By fitting the streaking curves (first moment of 
the final electron momentum distribution, white lines in Fig. la) to the analytic form of the 
IR vector potential Ai^{t + ts), we obtain the absolute delays shown in Fig. lb. All results 
are obtained by a nonlinear least-squares fit (error bars are present in all figures) of the first 
moment of the spectrogram taken along the laser polarization axis with an opening angle of 
10°. The spectrograms are obtained with a full numerical solution of the 3D TDSE without any 
approximations. The sign convention for ts ensures that positive values correspond to delayed 
emission, i.e. the electron 'feels' the vector potential at a later time. A negative time shift, or 
time advance, ts < 0, relative to the vector potential can be observed over a wide range of final 



electron energies (Fig. lb) corresponding to varying the XUV pulse frequency cjxuv- 

The question arises whether these time advances are related to the intrinsic quantum 
mechanical time shift corresponding to the scattering phases of the system. For this intrinsic 
time shift we use the single channel form of the Eisenbud-Wigner- Smith time shift tEws [12, 13] 
given by the energy derivative of the spectral phase 93 = arg{{po\z\il!i)), i.e. the group delay [19] 
of the electronic wave packet 

tEWS = (5) 

For photoionization in pure Coulomb potentials the phase ip is given by the Coulomb phase 
(Ti = aig\T{i+l—iZ /k)], where i is the angular momentum of the free electron, k its wavenumber, 
and Z is the charge of the remaining ion. However, the EWS time shift is only correctly defined 
for short-ranged potentials, so it comes as no great surprise that dai/dE does not agree with 
the time shift extracted from streaking (cf. Fig. lb). 

The physical interpretation of the extracted streaking time shifts requires a detailed inquiry 
into the assumptions underlying Eq. 4. The fundamental assumption of the streaking field 
as a probe is that the IR field does not distort the system under scrutiny. In the following 
we investigate the influence of the IR field on the entrance and exit channels. The original 
attosecond streaking model [3] is based on the strong-field approximation which assumes a 
'sudden' transition to the undisturbed momentum po caused by photoabsorption (Eq. 4) followed 
by a momentum shift by the IR field of the remaining streaking pulse. However, as the wave 
packet recedes from the ionic core, it propagates in the atomic (ionic) potential with a local 
momentum p{r) rather than the asymptotic momentum pQ. Thus, neglect of the deviation 
of the local from the asymptotic momentum is no longer valid for streaking in the Coulomb 
field (cf. [20]). Our classical trajectory Monte Carlo (CTMC) simulation (see [14] for details) 
gives streaking time shifts in almost perfect agreement with the full TDSE over a wide range 
of final-state energies for both H(ls) and He~''(ls) initial states (Fig. lb) and over two orders 
of magnitude in intensity (/ir ~ 10^*^ — lO^^W/cm^) of the streaking field (not shown). We 
therefore conclude that the energy dependence of the streaking time shifts results from classical 
Coulomb-laser coupling [14] (see also [20]). This has also been identified in one-dimensional 
studies using the eikonal approximation [21] and by an alternative interferometric method for 
extraction of time shifts and atomic phases ('RABITT') [15]. 

Initial state distortions by the IR field can be investigated by the constrained TDSE 
calculations (RIM) where the IR field explicitly only acts on the ionized part of the wave function. 
For energetically isolated initial states with small polarizability (H(ls), He~''(ls),...) the RIM 
results agree perfectly with the full TDSE results (see Fig. lb), therefore excluding initial- 
state polarization effects. However, if the initial state is not isolated, initial-state polarization 
can strongly affect the extracted delays. In Fig. 2 we show TDSE and RIM streaking results 
for low-lying states of He"'"(n£m) which are stable against tunnel ionization by the IR field 
(/jR < lO^^W/cm^). Within the full TDSE calculation we find remarkably strong initial-state 
dependent streaking field distortions with relative time shifts of up to 20 as between states with 
n = 2 and different quantum numbers^ £ and m. This pronounced streaking phase shift as a 
function of angular momentum differs from recent calculations [21, 22] which did not find a shift 
for states with (initially) vanishing dipole moment. However, for the simulations using the RIM 
the relative time shifts between Is and 2s disappear completely, and the time shifts for 2pQ and 
2pi almost agree (note that for initial p states the angular distributions are more complex and ts 

^ The strong variation of ts with angular momentum I is not specific to the degeneracy of the hydrogenic 
He^(n = 2) manifold. By choosing an atomic model potential that breaks the SO{4) symmetry we have verified 
that the streaking time shift persists for non-hydrogenic £ manifolds, thereby excluding degeneracy effects as the 
origin. Moreover, transient inter-shell coupling was found to dominate dynamical polarization over intra-shell 
mixing. 
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Figure 2: Streaking time shifts ts extracted from quantum mechanical streaking simulations 
with full TDSE (a) and RIM (b) for different initial states. The remaining difference between 
initial s and p states in the RIM results (b) can be attributed to the short-ranged scattering 
phase shift induced by the centrifugal potential (see Fig. 3). 
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Figure 3: Streaking time shifts t^iM. (lines) 
for excited initial states (2s and 2po) in 
He"*" (see Fig. 2b). The difference between 
tRiM(2po) and tRiM(2s) is given by A^ews = 
^ arg ( (po k 1 02po ) - M ( k I '^25 ) ) , which is 
given by the different centrifugal potential seen 
by the outgoing electron (V/ = /(/ + l)/(2r^) 
with I being the angular momentum of the free 
electron, / = p for initial s states and / = s or 
Z = d for initial p states) . 



becomes also more sensitive to the opening angle.) We therefore can rule out a state- dependent 
exit channel distortion of the IR field. 

The remaining difference between s and p states, after we have removed initial state 
distortion effects, can be finally attributed to the EWS time shift Atews = ^ arg((po|-z|02po) ~ 
^ arg((pokl</'2s)) as is shown in Fig. 3. Here ^ arg((po|-2|'/'2po) still includes the Coulomb 
phase shift o";, however, for different angular momenta I of the final state. The effective potential 
includes the long range Coulomb potential —Z/r, which is the same for both initial states, and 
the different centrifugal potential Vi = l{l + l)/(2r^). Thus we find that potentials are 
sufficiently short-ranged so that the corresponding scattering phases are accessible by streaking. 

We have shown in [14] that for pure short-ranged Yukawa potentials as well as for a short- 
ranged admixture to the Coulomb potential, the EWS-delay associated with the short-range part 
of the scattering phase agrees remarkably well with the observed streaking time shift. Along the 
same lines it is possible to describe streaking time shifts in static atomic model potentials by 
the Coulomb-laser coupling time shift of the asymptotic Coulomb field and the EWS time delay 
of the remaining short-ranged parts of the model potential as long as initial-state polarization 
can be neglected. 



4. Summary 

We have shown that Eisenbud-Wigner-Smith (EWS) time shifts (or energy variation of the scat- 
tering phase) for short-ranged potentials become accessible by attosecond streaking provided 
both initial-state dependent entrance channel and final-state exit channel distortions are prop- 
erly accounted for. For Coulomb potentials the coupling between the IR streaking field and the 
Coulomb field which depends on the final energy of the free electron dominates the extracted 
streaking time shift but can be accounted for classically. In addition we identified considerable 
state-dependent time shifts for easily polarizable initial states which are of quantum mechanical 
origin. Accounting for polarization of the initial state, the remaining difference of time delays 
between ionization from states with different angular momentum can be related to the EWS 
delay of the centrifugal potential. 
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